Quantum criticality and minimal conductivity in graphene with long-range disorder 
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We consider the conductivity a^x of graphene with negligible intervalley scattering at half filling. 
We derive the efi'ective field theory, which, for the case of a potential disorder, is a symplectic-class 
(j-model including a topological term with = n. As a consequence, the system is at a quantum 
critical point with a universal value of the conductivity of the order of e'^/h. When the effective 
time reversal symmetry is broken, the symmetry class becomes unitary, and a^x acquires the value 
characteristic for the quantum Hall transition. 



Recent breakthrough in graphene fabrication [I] and 
subsequent transport experiments [2] revealed remark- 
able electronic properties of this material. One of the 
most striking experimental observation is the minimal 
conductivity axx of order e^//i observed in undoped sam- 
ples and staying almost constant in a wide range of tem- 
peratures T from 300 K down to IK. This behavior 
should be contrasted to well-established results on the 
conductivity of two-dimensional (2D) systems where lo- 
calization effects drive the system into insulating state 
at low T 0, An apparently T-independent value of 
O'xx ~ e'^/h suggests that the system is close to a quan- 
tum critical point and calls for a theoretical explanation. 

One particular class of randomness when this scenario 
is realized, namely, the chiral disorder, was analyzed in 
detail in Ref. 5] (see also Ref. fi^]). It was shown that, 
if one of the chiral symmetries of clean graphene is pre- 
served by disorder, the conductivity at half filling is not 
affected by localization and is equal to (4/7r)e^//i (up 
to small corrections). While various types of random- 
ness in graphene (in particular, dislocations, ripples, or 
strong point-like defects) do belong to the chiral type, 
the experimentally observed value of axx is larger by a 
factor ^ 3, suggesting a different type of criticality. In 
this paper we consider another broad class of random- 
ness in graphene — long-range disorder. This case has 
a particular experimental relevance if the conductivity is 
dominated by charged impurities; the ripples 0, @] be- 
long to this class of randomness as well. Numerical sim- 
ulations of graphene with long-range random potential 
@j fiot provide an evidence in favor of a scale-invariant 
conductivity. 

The low-energy electron spectrum of clean graphene 
split into two degenerate valleys. The characteristic fea- 
ture of the long-range disorder is the absence of valley 
mixing due to the lack of scattering with large momen- 
tum transfer. This allows us to describe the system in 
terms of a single- valley Dirac Hamiltonian with disorder. 



H = wocrk + (TpV^(r). 



Here vq ~ 10* cm/s is the Fermi velocity. The four Pauli 
matrices cr^ (with ctq = 1) operate in the space of two- 
component spinors reflecting the sublattice structure of 



the honeycomb lattice, cr — {ai,a2}, and disorder in- 
cludes random scalar (Vq) and vector (Vi,2) potentials 
and random mass (V3). The Hamiltonian ([T]) was consid- 
ered in Ref. [ll| as a model for quantum Hall transition. 

To derive the field theory, we introduce a vector su- 
perfield ip with 2x2x2 = 8 components: the matrix 
structure of H in the sublattice space is complemented 
by the boson-fermion (BF) and the retarded-advanced 
(RA) structures. Assuming for simplicity Gaussian 5- 
correlated disorder distribution, we get the action 



, (2) 



with e = e + iOA where A — diag{l, — 1}^^- Assum- 
ing the isotropy, the disorder is described by three cou- 
plings ao = {V;i)/2TTvl, a± = {V^ + Vi)/2TTvl, and 
az = {V^)/2ttvq. On short (ballistic) scales the param- 
eters of di]) are renormalized 0, H, [ll|, [l3| ; the effective 
theory on longer scales is the non-linear sigma model [l3| . 

The clean single- valley Hamiltonian ([1]) obeys the ef- 
fective time-reversal (TR) invariance H = (J2H*(J2. This 
symmetry (denoted as Tl in Ref. 0) is not the physi- 
cal TR symmetry: the latter interchanges the nodes and 
is of no significance in the absence of inter-node scatter- 
ing. If the only disorder is random scalar potential, the 
TR invariance is not broken and the system falls into the 
symplectic symmetry class (AH) [3, 4|,|ll|. The standard 
realization of such symmetry is a system with spin-orbit 
coupling; in the present context the role of spin is played 
by the sublattice space. 

We start with a more generic case of the unitary sym- 
metry (class A) . The TR invariance is broken as soon as a 
(either random or non-random) mass or vector potential 
is included, in addition to the scalar potential. We find 
it instructive to present the derivation for a system with 
a non-zero mass term mfis. Decoupling the ip'^ term by 
a supermatrix field Q and integrating out ip, we get 



(1) S[Q] = Str 



dvrwgao 



ln( e — mcs -I- ivocrV + ijQ 



(3) 

where Str includes the matrix supertrace and the spatial 
integration, 7 = 1/2t, and r is the mean free time. 
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The saddle-point approximation 1^ reduces the set of 
Q to the conventional manifold of the unitary a-model; 
the relevant Q's are 4x4 supermatrices operating in RA 
and BF spaces and satisfying the constraint = 1. The 
low-energy modes describe slow spatial variation of Q on 
this manifold, and the effective theory is the result of 
the gradient expansion of the action ([3]) in these modes. 
As we show, this expansion is highly non-trivial due to 
anomalies in the theory of Dirac fermions |15l |. which 



induce a topological contribution to the cr-model !l0|. 

We first perform the gradient expansion of the real part 
of the action ([3]) 



Si[Q] = (l/2)Strln(G;iGli 



(4) 



Here the matrix Green functions defined as = e — 
m(73 -I- icr V ± Z7A are diagonal in RA space with retarded 
and advanced Green functions as their elements, G+ = 
diag{G^,G^}, G_ = diag{G^,G^}. Expanding Eq. Q 
to the second order in VQ and using the identity G^ — 
= —2i'fG^G^, we get the familiar gradient term. 



Si[Q] - (a,,/8)Str(VQ)' 



(5) 



The factor axx in Eq. ([4]) is identified as the dimcnsionless 
(in units /h) longitudinal conductivity given by 

cJx. - -(t^^/2)Tr[(G« - G^)a4G« - G^)a, 

= (l/27r) [1 + /(e, m){e' + 7' - m')/2^] , (6) 

where we introduced the notation 

f{x, y) = x^^{arctan[(.x + y)/j] + arctan[(a: — y)/'^]}- 

The calculation of the imaginary part iS2[Q] is much 
more subtle. We use the representation Q — T^^AT and 
cycle the matrices under the supertrace. The resulting 
expression depends on the vector u = TVT"^, 

iS2[Q] = (l/2)Str[ln(G;^+iwocru)-ln(Gl^+ii;ocru) 

The permutation of matrices leading to this formula is 
equivalent to a rotation of fermion fields, ip <—>■ Tip, in 
Eq. (12). This is not an innocent procedure in view of 
quantum anomaly 17|. However, such anomalous con- 
tributions from the two logarithms cancel in iS2[Q]- We 
proceed with expanding 15*2 [Q] in powers of u. The first 
two terms of this expansion are 

^4^)= (it,o/2)Str[A(G^-G^)cru] =7rStr(AJu), (7) 



is: 



(2) _ ^xy^afi 



2 Str(u„Au/3)= ^Str(gV,gVyQ). (8) 



The factors in Eqs. ([7]), ([5]) are the current spectral den- 
sity J(r) and the classical part of Hall conductivity [3| 



' xy 



{vll2) Tr[(G« + G^)<7x{G'' - G^)ay 

= -(m/27r)/(£,m). (9) 



The net current, and hence the linear term ([7]), is absent 
in the bulk of the system. It is incorrect, however, to drop 
this term. The contribution iS2^^ accounts for the edge 
current and gives the quantum part of the Hall conduc- 
tivity [lH. Prior to considering it, we have to establish 
boundary conditions (BC) for the Hamiltonian ([T]). 

Generically, BC in realistic graphene mix states from 
the two valleys of the spectrum. We can stay, how- 
ever, within the single-valley model and assume an in- 
finite mass M cx) at the boundary of the sample [l^l ■ 
Localization effects described by the cr-model occur in 
the bulk and hence are insensitive to particular BC. We 
thus assume that m(r) changes from a constant value m 
inside the sample to another, large value M outside it. 
The gradient of mass is not zero near the edge only. We 
further assume that the mass variation is slow on the 
scale of the electron mean free path but fast compared 
to (T-model length scales. This allows us to perform an 
expansion of the Green functions in Eq. ([7]) in Vm. With 
the help of identity [r, G] = ivoGcrG, we obtain 



Ja(r) = -ivl/2TT) Vpm tr a^G^asG'^apG'^ 

— <TaG^<T2,G^<TfjG^ 



The emerged trace is a mass derivative of the quantum 
part Gxy of Hall conductivity [3]; its direct calculation 
and subsequent integration with respect to m yields 

ag^~{m/2T:) f{m,e). (11) 

Substituting (llOp in ([7]) and integrating over the bound- 
ary strip, we express the term ([7]) as an integral along 
the edge and then apply the Stokes theorem: 



ea/3 



^iy{^) - ^iy{M)\ Str(AV„U0) 

sgn M ^ 



II 

4 \ 



Str(QV,gVyQ). (12) 



To derive the last expression, we have used the identity 
EQ/sVctW^ = eapupUa and the value of cr^^ in the limit 
of infinitely large M . The same result is obtained if one 
uses alternative BC introducing the second node with the 
large mass M. In that case, (l/2)sgnM will enter the 
action as the contribution of the second node to axy (cf. 
Refs. [il|,[l3]). Both contributions to 15*2 [Q], Eqs. ^ and 
(fT2)) . contain the functional Str(gVa;QVyQ) = 16z7r7V[Q] 
that is a well-known topological invariant on the cr-model 
manifold [Hi; its possible values are integer multiples of 
IBtti. The imaginary part of the action is defined up to 
a multiple of 27ri. Thus the sign of M is irrelevant, as 
expected: the bulk theory should not be sensitive to BC. 

Collecting all the contributions, we get the cr-model 
action for the single-node Dirac fermions with mass m: 



S[Q] = \ Str 



'xy 



(13) 
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The topological term is equal to i9N[Q], with the angle 



9 = 2Tiaxy 



In graphene the mass m is absent, so that 



<^xy — 0. Thus, the topological angle is, 9 — n. The the- 
ory (jl3p is then exactly on the critical line of the quantum 
Hall transition [l8| . in agreement with the arguments of 
Ref. Thus the graphene with a generic (TR break- 
ing) long-range disorder is driven into the quantum Hall 
critical point, with the conductivity 4a^ (the factor 4 
accounts for the spin and valley degeneracy). The value 
alj is known to be in the range ~ 0.5 — 0.6 from nu- 
merical simulations . A schematic scaling function in 
this case is shown in Fig. [1^. While formally this result 
holds for any energy e, in reality it only works near half 
filling (where the bare conductivity is ~ e^//i); for other 
e the quantum Hall critical point would only be reached 
for unrealistic temperatures and system sizes. 

Let us now turn to the case of preserved TR invariance, 
describing in particular charged impurities. The system 
belongs then to the symplectic symmetry class AH. The 
derivation of the cr-model starts with the doubling of %p 
variables accounting for the TR symmetry 13]. Then 
Q is 8 X 8 matrix obeying an additional constraint of 
charge conjugation Q = Q. The real part of the action is 
calculated in the same way as in the unitary case, yielding 
Eq. 11]) with an additional factor 1/2. 

Since the partition function of the symplectic model 
is real, the imaginary part of the action S'2[Q] can take 
one of the two possible values, or tt. The discrete- 
ness of 5*2 [Q] suggests that it again should be propor- 
tional to a topological invariant on the cr-model manifold. 
A non-trivial topology may arise only in the compact 
(fermion) sector of Q. The corresponding target space is 
M.F — 0{4n)/0{2n) x 0(2n), where n is the number of 
fermion species. While for the conventional ("minimal") 
tj- model n = 1, larger values will arise if one considers 
higher-order products of Green functions. The topologi- 
cal invariant takes values from the homotopy group [2l| 



TT2[Ml 



Z X Z, 



TT2[Mf\ 



(14) 



The homotopy group in the case n = 1 is richer than 
for n > 2. Nevertheless, S2[Q] may take only two non- 
equivalent values. Hence only a certain Z2 subgroup 
of the whole Z x Z comes into play as expected: the phase 
diagram of the theory should not depend on n. 

To demonstrate the emerging topology explicitly and 
to calculate the topological invariant, let us analyze the 
case n = 1 in more detail. The generators of the com- 
pact sector are Hermitian skew-symmetric 4x4 matrices 
anticommuting with A = p^: P1T2, P2T0, P2T1, and P2T3. 
Here pi and are Pauli matrices in RA and TR space 
respectively. These generators split into two mutually 
commuting pairs, each generating a 2-sphere ("diffuson" 
and "Cooperon" sphere). Simultaneous inversion of both 
spheres leaves Q intact. Hence the compact sector of the 
model is the manifold (5*2 x S'2)/Z2. Thus two topological 
invariants, iVi.2[Q], counting the covering of each sphere. 
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FIG. 1: Schematic scaling functions for (a) unitary and (b) 
symplectic universality class with topological term 9 — tt. 
The ordinary case S = is shown by dashed lines. 



emerge in accordance with Eq. (jl4|) . The most general 
topological term is 16*2 = i9iNi + i92N2. Due to the TR 
symmetry, the action is invariant under interchanging the 
diffuson and Cooperon spheres, which yields 9i = 92 = 9 
where 9 is either or tt. The explicit expression for n = 1 
topological term can be written using u = T\/T~^ 

^S2[Q] = ^Str[(A±l)r2«„U/3] = t7T{Ni[Q] + N2m, 

o 

yielding 9 = n. The sign ambiguity here does not af- 
fect any observables. If the TR symmetry is broken, the 
Cooperon modes are frozen and the manifold is reduced 
to a single diffuson sphere with iS2 [Q] — iirNi [Q] . 

The ordinary symplectic theory with no topological 
term exhibits a metal-insulator transition at a^p « 1.4 
[2^ . If the conductivity is smaller than this critical value 
the localization drives the system into insulating state, 
while in the metallic phase, a > <J*gp, antilocalization 
occurs. Using the analogy with the quantum Hall tran- 
sition in the unitary class, we argue that the topological 
term with 9 = tt suppresses localization effects when the 
conductivity is small, leading to appearance of a new 
attracting fixed point at cr^* . The position of the metal- 
insulator transition, a^p, is also affected by the topolog- 
ical term. However, we believe that its change is neg- 
ligible: the instanton correction to the scaling function 
at large conductivity is exponentially small 18], and the 
value of the exponential factor e^*'^'^ is still extremely 
small at tJ = cr^p- A plausible scaling of the conductivity 
in the symplectic case with 6* = tt is sketched in Fig. [T] 
The existence and position of the new critical point can 
be verified numerically. Recent simulations of graphene 

[13] indeed demonstrate the stability of the conductiv- 
ity in the presence of long-range disorder. Of course, in 
reality there will be always a weak inter-valley scatter- 
ing, which will establish the localization and lowest T, in 
agreement with Q. However, the approximate quantum 
criticality will hold in a parametrically broad range of T. 

Finally, we discuss a connection between our findings 
and recent results on the quantum spin Hall (QSH) eS'ect 
in systems of Dirac fermions with spin-orbit coupling [2^ ] , 
which in the presence of random potential also belong to 
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the symplectic symmetry class. Such systems were found 
to possess two distinct insulating phases, both having a 
gap in the bulk electron spectrum but differing by the 
edge properties. While the normal insulating phase has 
no edge states, the spin-Hall insulator is characterized 
by a pair of mutually time-reversed spin-polarized edge 
states penetrating the bulk gap. The transition between 
these two phases is driven by Rashba spin-orbit coupling 
strength. The existence of the edge states was attributed 
to certain Z2 topological index different from one stud- 
ied above. This topological order is robust with respect 
to disorder, even if the latter mixes the valleys. The 2D 
cr-model is insensitive to the edge properties of the sam- 
ple and does not capture the difference between the two 
insulating phases. A suitable effective theory is the one- 
dimensional (ID) (T-model for the edge states. The cor- 
responding topological index characterizes the homotopy 
group tti(A4f)- This group is again Z2 that enables a 
0-term with 9 equal to or tt. The topological term with 
9 = TT is present if the number of channels is odd. Then 
the conductivity of the ID system equals e'^/h in the 
long- length limit psj : one conducting channel survives 
localization. This is what happens in QSH systems when 
a pair of edge states is not localized [2J]. In the pres- 
ence of disorder, such systems will possess three phases: 
metal, normal insulator, and QSH insulator. We expect 
that generically there will be a transition between the 
latter two. The critical theory discussed in the present 
work (2D symplectic cr-model with 9 — n) should then 
describe this QSH transition. 

In summary, graphcnc with long-range disorder shows 
quantum criticality at half filling. If the effective TR 
symmetry of the single-valley system is preserved (e.g. 
when Coulomb scatterers are the dominant disorder), the 
relevant theory is the symplectic tx-model with topologi- 
cal angle 9 = n and the minimal conductivity takes a uni- 
versal value 4(7 J* . If the TR symmetry is broken (e.g. by 
effective random magnetic field due to ripples), the sys- 
tem falls into the universality class of the quantum Hall 
critical point, with another universal value axx = 4(t^. 
We have argued that the symplectic critical point de- 
scribes also the QSH transition [j^ ]. 
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